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Abstract 

Correlations between the elliptic or triangular flow coefficients v m (m = 2 or 3) and other flow 
harmonics v n (n = 2 to 5) are measured using v / s ~ = 2.76 TeV Pb+Pb collision data collected in 
2010 by the ATLAS experiment at the LHC, corresponding to an integrated lumonisity of 7 /ub -1 . The 
v m -v n correlations are measured in midrapidity as a function of centrality, and, for events within the 
same centrality interval, as a function of event ellipticity or triangularity defined in a forward rapidity 
region. For events within the same centrality interval, v 3 is found to be anticorrelated with v 2 and this 
anticorrelation is consistent with similar anticorrelations between the corresponding eccentricities e 2 
and e 3 . On the other hand, it is observed that v 4 increases strongly with v 2 , and v 5 increases strongly 
with both v 2 and v 3 . The trend and strength of the v m -v n correlations for n = 4 and 5 are found to 
disagree with e m -e n correlations predicted by initial-geometry models. Instead, these correlations are 
found to be consistent with the combined effects of a linear contribution to v n and a nonlinear term 
that is a function of v\ or of v 2 v 3 , as predicted by hydrodynamic models. A simple two-component fit is 
used to separate these two contributions. The extracted linear and nonlinear contributions to v 4 and 
v 5 are found to be consistent with previously measured event-plane correlations. 
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Measurement of the correlation between flow harmonics of different order in lead—lead 
collisions at y/s NN = 2.76 TeV with the ATLAS detector 

The ATLAS Collaboration 

Correlations between the elliptic or triangular flow coefficients v m {jn = 2 or 3) and other flow 
harmonics v n (n = 2 to 5 ) are measured using y/s NN = 2.76 TeV Pb+Pb collision data collected in 
2010 by the ATLAS experiment at the LHC, corresponding to an integrated lumonisity of 7 pb _1 . 

The v m -v n correlations are measured in midrapidity as a function of centrality, and, for events within 
the same centrality interval, as a function of event ellipticity or triangularity defined in a forward 
rapidity region. For events within the same centrality interval, v$ is found to be anticorrelated with 
V 2 and this anticorrelation is consistent with similar anticorrelations between the corresponding 
eccentricities £2 and £ 3 . On the other hand, it is observed that 14 increases strongly with v 2 , and v$ 
increases strongly with both V2 and 113 . The trend and strength of the Vm~v„ correlations for n = 4 
and 5 are found to disagree with e m -£n correlations predicted by initial-geometry models. Instead, 
these correlations are found to be consistent with the combined effects of a linear contribution to 
v n and a nonlinear term that is a function of d 2 or of V2V3 , as predicted by hydrodynamic models. 

A simple two-component fit is used to separate these two contributions. The extracted linear and 
nonlinear contributions to va and v$ are found to be consistent with previously measured event-plane 
correlations. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 


Heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) create 
hot and dense matter that is thought to be composed of strongly coupled quarks and gluons. The distribution 
of this matter in the transverse plane is both non-uniform in density and asymmetric in shape [1, 2]. The matter 
expands under large pressure gradients, which transfer the inhomogenous initial condition into azimuthal anisotropy of 
produced particles in momentum space [3, 4]. Hydrodynamic models are used to understand the space-time evolution 
of the matter by comparing predictions with the measured azimuthal anisotropy [5-7] . The success of these models in 
describing the anisotropy of particle production in heavy-ion collisions at RHIC and the LHC [8-14] places significant 
constraints on the transport properties (such as the ratio of shear viscosity to entropy density) and initial conditions 
of the produced matter [15-20]. 

The azimuthal anisotropy of the particle production in each event can be characterized by a Fourier expansion of 
the corresponding probability distribution P((f>) in azimuthal angle <f) [3, 21 ], 
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where v n and 4>„ are the magnitude and phase (also known as the event plane or EP) of the n th -order harmonic 
flow, and is real by construction. The presence of harmonic flow has been related to various moments of shape 
configurations of the initially produced fireball. These moments are described by the eccentricity vector e n calculated 
from the transverse positions (r, <fi) of the participating nucleons relative to their center of mass [4, 16]: 
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where (...) denotes an average over the transverse position of all participating nucleons, and e„ and „ (also known 
as the participant plane or PP) represent the magnitude and orientation of the eccentricity vector, respectively. 
The eccentricity vectors characterize the spatial anisotropy of the initially produced fireball, which drives the flow 
harmonics in the final state. 

According to hydrodynamic model calculations, elliptic flow V 2 and triangular flow V 3 are the dominant harmonics, 
and they are driven mainly by the ellipticity vector e 2 and triangularity vector e 3 of the initially produced fireball [ 22 , 
23]: 
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This proportionality is often quantified by a ratio 

k n = v n /e n , n = 2 or 3, (4) 


where the linear response coefficients k n are found to be independent of the magnitude of e„ but change with central¬ 
ity [22, 24], 

The origin of higher-order (n > 3) harmonics is more complicated; they arise from both e n and nonlinear mixing of 
lower-order harmonics [20, 23, 25]. For example, an analytical calculation shows that the V 4 signal comprises a term 
proportional to £4 (linear response term) and a leading nonlinear term that is proportional to [23, 26]: 

u 4 e i4 * 4 = ao^+fl! (e ^ 2 * 2 ) 2 + ... 

= c 0 e l44>i + ci (v 2 e l24>2 ) 2 + ... , (5) 


where the second line of the equation follows from Eq. (3), and Co = aoE 4 denotes the linear component of u 4 and 
coefficients a 0 , a\ and C\ are weak functions of centrality. The nonlinear contribution from v 2 is responsible for the 
strong centrality dependence of the correlation between $2 and $ 4 observed by the ATLAS Collaboration [14] in 
Pb+Pb collisions. In a similar manner, the v§ signal comprises a linear component proportional to £5 and a leading 
nonlinear term involving v 2 and v 3 [23, 26]: 
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This decomposition of the v 3 signal explains the measured EP correlation involving 4> 2 , *I > 3 and $5 [14]. 

Due to fluctuations of nucleon positions in the initial state, e n and v n vary from event to event, which can be 
described by probability distributions p(e n ) and p{v n ). Recent measurements by the ATLAS Collaboration [13] show 
that the distributions p(v n ) are very broad: even for events in a very narrow centrality interval, v 2 and v 3 can fluctuate 
from zero to several times their mean values. If events with different v 2 or v 3 values could be selected cleanly, one 
would be able to control directly the relative sizes of the linear and nonlinear contributions to V 4 and v 3 in Eqs. (5) 
and ( 6 ), and hence provide an independent method of separating these two contributions. Such an event-shape 
selection method has been proposed in Refs. [27, 28], where events in a narrow centrality interval are further classified 
according to the observed ellipticity or triangularity in a forward rapidity region. These quantities are estimated 
from the “flow-vector” q m (to = 2 and 3) as described in Section IV A. This classification gives events with similar 
multiplicity but with very different ellipticity or triangularity. By measuring the v n and v m in a different rapidity 
window for each q m event class, the differential correlation between v m and v n can be obtained in an unbiased way 
for each centrality interval, which allows the separation of the linear and nonlinear components in t ) 4 and v 3 . The 
extracted linear component of V 4 and v 3 can then be used to understand the collective response of the medium to the 
initial eccentricity of the same order, using an approach similar to Eq. (4). 

In addition to separating the linear and nonlinear effects, the correlation between v m and v n is also sensitive to any 
differential correlation between e m and e n in the initial state. One example is the strong anticorrelation between e 2 
and £ 3 predicted by the Monte Carlo (MC) Glauber model [28, 29]. A recent transport-model calculation shows that 
this correlation survives the collective expansion and appears as a similar anticorrelation between v 2 and v 3 [28]. 

In this paper, the correlations between two flow harmonics of different order are studied using the event-shape 
selection method. The ellipticity or triangularity of the events is selected based on the q 2 or q 3 signal in the forward 
pseudorapidity range of 3.3 < |? 7 | < 4.8 . 1 The values of v n for n = 2 to 5 are then measured at midrapidity |? 7 | < 2.5 
using a two-particle correlation method, and the correlations between two flow harmonics are obtained. The procedure 
for obtaining v n in this analysis is identical to that used in a previous ATLAS publication [ 11 ], which is also based on 
the same dataset. The main difference is that, in this analysis, the events are classified both by their centrality and 
by the observed q 2 or q 3 at forward pseudorapidity. Most systematic uncertainties are common to the two analyses. 


1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the detector and 
the 2 -axis along the beam pipe. The ir-axis points from the IP to the center of the LHC ring, and the ;/-axis points upward. Cylindrical 
coordinates (r, 0 ) are used in the transverse plane, cf> being the azimuthal angle around the beam pipe. The pseudorapidity is defined in 
terms of the polar angle 6 as 77 = — lntan( 0 / 2 ). 
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II. ATLAS DETECTOR AND TRIGGER 

The ATLAS detector [30] provides nearly full solid-angle coverage of the collision point with tracking detectors, 
calorimeters and muon chambers. All of these are well suited for measurements of azimuthal anisotropies over a 
large pseudorapidity range. This analysis primarily uses two subsystems: the inner detector (ID) and the forward 
calorimeter (FCal). The ID is contained within the 2 T field of a superconducting solenoid magnet and measures 
the trajectories of charged particles in the pseudorapidity range \r/\ < 2.5 and over the full azimuth. A charged 
particle passing through the ID traverses typically three modules of the silicon pixel detector (Pixel), four double¬ 
sided silicon strip modules of the semiconductor tracker (SCT), and a transition radiation tracker for |? 7 | < 2. The 
FCal consists of three sampling layers, longitudinal in shower depth, and covers 3.2 < |?y| < 4.9. The energies in 
the FCal are reconstructed and grouped into towers with segmentation in pseudorapidity and azimuthal angle of 
Ar] x A <f> ss 0.2 x 0.2. In heavy-ion collisions, the FCal is used mainly to measure the event centrality and event 
planes [11, 31]. In this analysis it is also used to classify the events in terms of q 2 or q% in the forward rapidity region. 

The minimum-bias trigger used for this analysis requires signals in two zero-degree calorimeters (ZDC) or either of 
the two minimum-bias trigger scintillator (MBTS) counters. The ZDCs are positioned at ±140 m from the collision 
point, detecting neutrons and photons with |? 7 | > 8.3, and the MBTS covers 2.1 < |? 7 | < 3.9 on each side of the nominal 
interaction point. The ZDC trigger thresholds on each side are set below the peak corresponding to a single neutron. 
A timing requirement based on signals from each side of the MBTS is imposed to remove beam backgrounds. 


III. EVENT AND TRACK SELECTION 

This analysis is based on approximately 7 /ib _1 of Pb±Pb data collected in 2010 at the LHC with a nucleon- 
nucleon center-of-mass energy s NN = 2.76 TeV. The offline event selection requires a reconstructed vertex and a 
time difference |Af| < 3 ns between signals in the MBTS trigger counters on either side of the interaction point 
to suppress noncollision backgrounds. A coincidence between the ZDCs at forward and backward pseudorapidity is 
required to reject a variety of background processes, while maintaining high efficiency for inelastic processes. Events 
satisfying these conditions are further required to have a reconstructed primary vertex with |z v t x | < 150 mm from the 
nominal center of the ATLAS detector. About 48 million events pass the requirements. 

The Pb±Pb event centrality [32] is characterized using the total transverse energy (SAt) deposited in the FCal 
over the pseudorapidity range 3.2 < |7y| < 4.9 at the electromagnetic energy scale [33]. From an analysis of this 
distribution after all trigger and event-selection requirements, the fraction of the inelastic cross-section sampled is 
estimated to be 98 ±2%. The uncertainty associated with the centrality definition is evaluated by varying the effect of 
trigger and event selection inefficiencies as well as background rejection requirements in the most peripheral FCal EAp 
interval [32]. The FCal EAt distribution is divided into a set of 5% percentile bins. A centrality interval refers to a 
percentile range, starting at 0% relative to the most central collisions. Thus the 0-5% centrality interval corresponds 
to the most central 5% of the events. An MC Glauber analysis [32, 34] is used to estimate the average number of 
participating nucleons, A part , for each centrality interval. These are summarized in Table I. Following the convention 
of heavy-ion analyses, the centrality dependence of the results in this paper is presented as a function of A part . 


TABLE I: The list of centrality intervals and associated values of the average number of participating nucleons Apart used in 
this analysis. The systematic uncertainties are taken from Ref. [32]. 

Centrality 0-5% 5-10% 10-15% 15-20% 20-25% 25-30% 30-35% 

A part 382 ± 2 330 ±3 282 ± 4 240 ± 4 203 ± 4 170 ± 4 142 ±4 
Centrality 35-40% 40-45% 45-50% 50-55% 55-60% 60-65% 65-70% 

Ap art 117 ±4 95 ± 4 76 ±4 60 ± 3 46 ± 3 35 ± 3 25 ± 2 


The harmonic flow coefficients v n are measured using tracks in the ID that are required to have transverse momentum 
Pt > 0.5 GeV and |? 7 | < 2.5. At least nine hits in the silicon detectors are required for each track, with no missing 
Pixel hits and not more than one missing SCT hit, taking into account the effects of known dead modules. In addition, 
the point of closest approach of the track is required to be within 1 mm of the primary vertex in both the transverse 
and longitudinal directions [31]. The efficiency e(pTj??) of the track reconstruction and track selection requirements 
is evaluated using simulated Pb±Pb events produced with the HIJING event generator (version 1.38b) [35]. The 
generated particles in each event are rotated in azimuthal angle according to the procedure described in Ref. [36] to 
give harmonic flow consistent with previous ATLAS measurements [11, 31]. The response of the detector is simulated 
using GEANT4 [37, 38] and the resulting events are reconstructed with the same algorithms that are applied to the 
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data. The absolute efficiency increases with px by 7% between 0.5 GeV and 0.8 GeV and varies only weakly for 
Px > 0.8 GeV. However, the efficiency varies more strongly with 77 and event multiplicity [31]. For p T > 0.8 GeV, it 
ranges from 72% at 77 « 0 to 57% for |? 7 | > 2 in peripheral collisions, while it ranges from 72% at 77 « 0 to about 42% 
for 177 ! > 2 in central collisions. 


IV. DATA ANALYSIS 


A. Event-shape selection 


The ellipticity or triangularity in each event is characterized by the so-called “flow vector” calculated from the 
transverse energy (Ex) deposited in the FCal [14, 39]: 
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where the weight Wj is the Et of the j th tower at azimuthal angle <f>j in the FCal. Subtraction of the event-averaged 
centroid (q m ) evts in Eq. (7) removes biases due to detector effects [40]. The angles are the observed event 
planes, which fluctuate around the true event planes <F m due to the finite number of particles in an event. A standard 
technique [41] is used to remove the small residual nonuniformities in the distribution of 4/^ s . These procedures are 
identical to those used in several previous flow analyses [11, 13, 14, 40]. To reduce the detector nonuniformities at 
the edge of the FCal, only the FCal towers whose centroids fall within the interval 3.3 < |? 7 | < 4.8 are used. 

The q m defined above is insensitive to the energy scale in the calorimeter. In the limit of infinite multiplicity it 
approaches the Ex-weighted single-particle flow: 
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Hence the q m distribution is expected to follow closely the v m distribution, except that it is smeared due to the 
finite number of particles. Figure 1 shows the distributions of (72 and q 3 in the 0-1% most central collisions. These 
events are first divided into ten q m intervals with equal number of events. Since the intervals at the highest and 
lowest q m values cover much broader ranges, they are further divided into five and two smaller intervals, respectively, 
resulting in a total of fifteen q m intervals containing certain fractions of events. Starting at the low end of the q m 
distribution, there are two intervals containing a fraction 0.05 (labeled 0.95-1 and 0.9-0.95), eight intervals containing 
0.1, three containing 0.025, one containing 0.015, and one containing 0.01 (this last interval spans the highest values 
of q m ). These fifteen intervals are defined separately for each 1% centrality interval, and are then grouped together 
to form wider centrality intervals used in this analysis (see Table I). For example, the first q m interval for the 0-5% 
centrality interval is the sum of the first q m interval in the five centrality intervals, 0-1%, 1-2%,..., 4-5%. The default 
analysis uses fifteen nonoverlapping q m intervals defined in Fig. 1. For better statistical precision, sometimes they are 
regrouped into wider q m intervals. 


B. Two-particle correlations 

The two-particle correlation analysis closely follows a previous ATLAS publication [11] where it is described in detail, 
so the analysis is only briefly summarized here. For a given event class, the two-particle correlation is measured as a 
function of relative azimuthal angle A = (f> a — 4>b and relative pseudorapidity A 77 = 773 , — 77 b ■ The labels a and b denote 
the two particles in the pair, which may be selected from different px intervals. The two-particle correlation function 
is constructed as the ratio of distributions for same-event pairs (or foreground pairs S(A(j) 1 A 77 )) and mixed-event 
pairs (or background pairs B(A(f>, A 77 )): 


C(A<j), A 77 ) 


SjA^A 77 ) 
B(A<j), A 77 )' 


(9) 


The mixed-event pair distribution is constructed from track pairs from two separate events with similar centrality and 
z v tx, such that it properly accounts for detector inefficiencies and nonuniformity, but contains no physical correlations. 
Charged particles measured by the ID with a pair acceptance extending up to | A 77 I = 5 are used for constructing the 
correlation function. 

This analysis focuses mainly on the shape of the correlation function in A <j>. A set of 1-D A <j> correlation functions 
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FIG. 1: (Color online) The distributions of the magnitude of the flow vector, q 3 (left panel) and q 3 (right panel), calculated 
in the FCal via Eq. (7) in the 1% most central collisions. The vertical lines indicate the boundaries of the fifteen q m ranges, 
each containing a fraction of events as indicated. 




FIG. 2: (Color online) The correlation functions C{A(j>) for pairs with |A»j| > 2 and 0.5 < pt < 2 GeV in 0-5% centrality. 
The correlation functions for events with the largest 10% and smallest 10% g m values are also shown for m = 2 (left panel) and 
m = 3 (right panel). The statistical uncertainties are smaller than the symbols. 


is built from the ratio of the foreground distributions to the background distributions, both projected onto A <f>: 


C( A0) 


/ S(A<p, Ai])dAi] 
f B(A(j), Ar/)dAij 


( 10 ) 


The normalization is fixed by scaling the number of the mixed-event pairs to be the same as the number of same-event 
pairs for 2 < |A^| < 5, which is then applied to all A 77 slices. 

Figure 2 shows the 1-D correlation functions for 2 < \Arj\ < 5 calculated in the low-px region (0.5 < p^ h < 2 GeV) 
in the 0-5% most central collisions. The correlation functions are also shown for events selected with the largest and 
smallest <72 values (left panel) or q 3 values (right panel). The magnitude of the modulation correlates strongly with the 
q m value reflecting the fact that the global ellipticity or triangularity can be selected by q 2 or q 3 in the forward rapidity 
interval. The correlation function for events with smallest q 2 or largest q 3 values shows a double-peak structure on 
the away-side (A <f> ~ n). This structure reflects the dominant contribution of the triangular flow under these q m 
selections. Similar double-peak structures are also observed in ultracentral Pb+Pb collisions without event-shape 
selection [11, 42]. 
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The 1-D correlation function in is then expressed as a Fourier series: 

C{ A(j>) = ^ ^1 + 2 ^2 v n,n cos (nAcp)^ . (11) 

The Fourier coefficients are calculated directly from the correlation function as v n ^ n = (cos (nA</>)). The single-particle 
azimuthal anisotropy coefficients v n are obtained via the factorization relation commonly used for collective flow in 
heavy ion collisions [11, 12, 43, 44]: 


v n,n{PT’Pr) = 


( 12 ) 


From Eq. (12), v n is calculated as: 


v n (pr) 


i,n(PT j Pt) /\J v n,n(PT>PT) 


(13) 


where p^ is simply denoted by pt from now on, and the default transverse momentum range for p\ is chosen to be 
0.5 < p\ < 2 GeV, where the hydrodynamic viscous corrections are not too large. The v n values obtained using 
this method measure, in effect, the root-mean-square (RMS) values of the event-by-event v n [43]. A detailed test 
of the factorization behavior was carried out [11, 12] by comparing the u n (px) obtained for different p\ ranges, and 
factorization was found to hold to within 10% for p\ < A GeV for the centrality ranges studied in this paper. 


C. Systematic uncertainties 

Other than the classification of events according to q m (m = 2 or 3), the analysis procedure is nearly identical to 
the previous ATLAS measurement [11] based on the same dataset. Most systematic uncertainties are the same, and 
they are summarized here. 

The correlation function relies on the pair acceptance function to reproduce and cancel the detector acceptance 
effects in the foreground distribution. A natural way of quantifying the influence of detector effects on and 
v n is to express the single-particle and pair acceptance functions as Fourier series [as in Eq. (11)], and measure the 
coefficients and v be *. The resulting coefficients for pair acceptance, v be ‘, are the product of two single-particle 
acceptances u bet,a and u bet,b . In general, the pair acceptance function in A <j> is quite flat: the maximum variation 
from its average is observed to be less than 0.001, and the corresponding values are found to be less than 

1.5 x 10~ 4 . These effects are expected to cancel to a large extent in the correlation function, and only a small 
fraction contributes to the uncertainties in the pair acceptance function. Three possible residual effects for v are 
studied in Ref. [11]: 1) the time dependence of the pair acceptance, 2) the effect of imperfect centrality matching, and 
3) the effect of imperfect z vtx matching. In each case, the residual values are evaluated by a Fourier expansion of 
the ratio of the pair acceptances before and after the variation. The systematic uncertainty of the pair acceptance is 
the sum in quadrature of these three estimates, which is Sv„ tn < 5 x 10 -6 for 2 < |A??| < 5. This absolute uncertainty 
is propagated to the uncertainty in v n , and it is the dominant uncertainty when v n is small, e.g., for v$ in central 
collisions. This uncertainty is found to be uncorrelated with the q m selection, and hence it is assumed not to cancel 
between different q m intervals. 

A further type of systematic uncertainty includes the sensitivity of the analysis to track selection requirements and 
track reconstruction efficiency, variation of v n between different running periods, and trigger and event selection. The 
effect of the track reconstruction efficiency was evaluated in Ref. [13]; the other effects were evaluated in Ref. [11]. 
Most systematic uncertainties cancel in the correlation function when dividing the foreground distribution by the 
background distribution. The estimated residual effects are summarized in Table II. Most of these uncertainties are 
expected to be correlated between different q m intervals. 

Finally, due to the anisotropy of particle emission, the detector occupancy is expected to be larger in the direction 
of the event plane, where the particle density is larger. Any occupancy effects depending on azimuthal angle may 
lead to a small angle-dependent efficiency variation, which may slightly reduce the measured v n coefficients. The 
magnitude of such an occupancy-dependent variation in tracking efficiency is evaluated using the HIJING simulation 
with flow imposed on the generated particles [13]. The reconstructed v n values are compared to the generated v n 
signal. The differences are taken as an estimate of the systematic uncertainties. These differences are found to be a 
few percent or less, and are included in Table II. Since this effect is proportional to the flow signal, it is expected to 
partially cancel between different q m ranges. 
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TABLE II: Relative systematic uncertainties on the measured v n due to track selection requirements, track reconstruction 
efficiency, variation between different running periods, trigger selection, consistency between true and reconstructed v n in 
HIJING simulation, and the quadrature sum of individual terms. Most of these uncertainties are correlated between different 
ranges of q m (m = 2 or 3). 



V 2 

V 3 

zq 

V 5 

( 7 m -dependent 

Track selection [%] 

0.3 

0.3 

1.0 

2.0 

yes 

Track reconstruction efficiency [%] 

0 . 1 - 1.0 

0.2 1.5 

0 .2-2.0 

D.3-2.5 

yes 

Running periods [%] 

0 .3-1.0 

0.7-2.1 

1.2-3.1 

2.3 

no 

Trigger [%] 

0 .5-1.0 

0.5-1.0 

0.5-1 

1.0 

yes 

MC closure and occupancy effects [%] 

1.0 

1.5 

2.0 

3.5 

yes 

Sum of above [%] 

1 .2-2.0 

1.8-3.2 

2.6-4.4 

4.7-5.4 



V. RESULTS 

A. Fourier coefficients v n and their correlations with q m 

Figure 3 shows the zq(px) for n = 2 to 5 extracted via Eq. (13) for events in the 20-30% centrality interval. The 
results show nontrivial correlations with both the q 2 (left column) or (right column) selections. In the case of the 
q 2 selection, the V 2 values are largest for events selected with the largest q 2l and smallest for events selected with the 
smallest ( 72 , with a total change of more than a factor of two. A similar dependence on <72 is also seen for zq(px) and 
zq (pt ) (two bottom panels). In contrast, the extracted Z 3 (px) values are anticorrelated with < 72 ; the overall change 
in V 3 (pt) is also significantly smaller (< 20% across the (72 range). In the case of the (73 selection, a strong positive 
correlation is observed for zq and zq, and a weak anticorrelation is observed for V 2 and zq. All these correlations are 
observed to be nearly independent of px, suggesting that the response of v n to the change in the event shape is largely 
independent of px- As a consistency check, the inclusive results without q m selection are compared with previously 
published results from Ref. [11]: the differences are less than 0.6% for V 2 and increase to 2-3% for higher harmonics, 
which are well within the systematic uncertainties quoted in Table II. 

Figure 4 shows the correlation between v n and q m for m = 2 (left column) and m = 3 (right column) in several 
centrality intervals in a low px range (0.5 < px < 2 GeV). Since the v n -q m correlation depends only weakly on px, this 
plot captures the essential features of the correlation between v n and q m shown in Fig. 3. Due to the finite number of 
particles in an event, the measured q m values fluctuate relative to the true values, diluting the correlations with v n . 
The influence of smearing on the <72 is much smaller than that for the 53 simply because the V 2 signal is much bigger 
than the V 3 signal. But since both the v m -q m and v n -q m correlations are measured, the results are presented directly 
as v m -v n correlations for various q m selections. The level of detail contained in the v m -v n correlation is controlled by 
the dynamic range of v m when varying the q m selection. This dynamic range depends strongly on event centrality. 
For example, in the 10-15% centrality interval, V 2 is varied by a factor of 3.1 by selecting on (72 and V 3 is varied by 
a factor of 2.4 by selecting on 53 . In the 40-45% centrality interval, however, due to stronger statistical smearing of 
q m , the zq and V 3 are only varied by a factor of 2.7 and 1.7, respectively. Hence the event-shape selection is precise 
in central and midcentral collisions and is expected to be less precise in peripheral collisions. 

In general, correlations v m -q m and v n -q m can be measured in different px ranges, and the derived v m -v n correlation 
can be categorized into three types: 1) the correlation between v m in two different px ranges, v m {p^}-v m {p\}, 2) 
the correlation between v m and another flow harmonic of different order v„ in the same px range, z m {px}-z n {px}, 
and 3) the correlation between v m and v n in different px ranges, v m {p^}-v n {p\.}. However, the Vm{p^}-v n {p^} 
correlation can be obtained by combining two correlations u m {pT}-Wm{Px} and v m {p^}-v n {p\.}, so it does not carry 
independent information. This paper therefore focuses on the first two types of correlation. 

The results for v m -v n correlations are organized as follows. Section VB presents correlations of V 2 or V 3 between two 
different px ranges. The V 2 -V 3 correlations are discussed in Sec. V C. This is followed by tq-'tq and Z3-Z4 correlations 
in Sec. VD and zq-zq and zq-zq correlations in Sec. VE, where a detailed analysis is performed to separate the linear 
and nonlinear components of zq and zq. The eccentricity scaling behavior of the extracted linear component of v n is 
presented in Sec. VF. 
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FIG. 3: (Color online) The harmonic flow coefficients v u (pt) in the 20-30% centrality interval for events selected on either 
q 2 (left column) or 53 (right column) for n = 2 (top row), n = 3 (second row), n = 4 (third row) and n = 5 (bottom row). 
They are calculated for reference px of 0.5 < Px < 2 GeV [Eq. (13)]. The top part of each panel shows the u n (px) for events in 
the 0-0.1, 0.1-0.2, 0.7-0 .8 and 0.9-1 fractional ranges of q m (open symbols) as well as for inclusive events without q m selection 
(solid symbols). The bottom part of each panel shows the ratios of the v u (pt) for g m -selected events to those obtained for all 
events. Only statistical uncertainties are shown. 
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FIG. 4: (Color online) The correlations between v n and 52 (left column) and <73 (right column) in four centrality intervals with 
n = 2 (top row), n = 3 (second row), n = 4 (third row) and n = 5 (bottom row), where v n is calculated in 0.5 < pt < 2 GeV. 
Only statistical uncertainties are shown. The lines connecting data points are for guidance only. 
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FIG. 5: (Color online) The correlation of the v m between 0.5 < pt < 2 GeV (x-axis) and 3 < pt < 4 GeV (j/-axis) for m = 2 
(left panel) and m = 3 (right panel). The v m values are calculated for fourteen 5% centrality intervals in the centrality range 
0-70% without event-shape selection. The data points are connected to show the boomerang trend from central to peripheral 
collisions as indicated. The error bars and shaded boxes represent the statistical and systematic uncertainties, respectively. 
These uncertainties are often smaller than the symbol size. 


B. Correlation of V2 or V3 between two different px ranges 

Figure 5 shows the correlation of v m for m = 2 (left panel) or m = 3 (right panel) between two px ranges for 
various centrality intervals. The z-axis represents v m values in the 0.5 < px < 2 GeV range, while the p-axis 
represents v m values from a higher range of 3 < px <4 GeV. Each data point corresponds to a 5% centrality 
interval within the overall centrality range of 0-70%. Going from central collisions (left end of the data points) to the 
peripheral collisions (right end of the data points), v m first increases and then decreases along both axes, reflecting 
the characteristic centrality dependence of v m , well known from previous flow analyses [10, 11]. The rate of decrease is 
larger at higher px, resulting in a “boomerang-like” structure in the correlation. The stronger centrality dependence 
of v m at higher px is consistent with larger viscous-damping effects expected from hydrodynamic calculations [45] . 

In the next step, events in each centrality interval are further divided into q m intervals, as described in Sec. IV A. 
With this further subdivision, each data point in Fig. 5 turns into a group of data points, which may follow a 
different correlation pattern. These data points are shown in Fig. 6 (markers) overlaid with the overall centrality 
dependence prior to the event-shape selection from Fig. 5 (the “boomerang”). For clarity, the results are shown only 
for seven selected centrality intervals. Unlike the centrality dependence, the v m correlation within a given centrality 
interval approximately follows a straight line passing very close to the origin. The small nonzero intercepts can be 
attributed to a residual centrality dependence of the v m -v m correlation within the finite centrality intervals used. This 
approximately linear correlation suggests that, once the event centrality or the overall event multiplicity is fixed, the 
viscous-damping effects on v rn change very little with the variation of the event shape via q m selection. The influence 
of viscosity on flow harmonics is mainly controlled by the event centrality (or the overall system size). 


C. V2—V3 correlation 

Figure 7(a) shows the centrality dependence of the correlation between V 2 and V 3 measured in 0.5 < px < 2 GeV. 
The boomerang-like structure in this case reflects mostly the fact that V 3 has a much weaker centrality dependence 
than V 2 [11]. Figure 7(b) overlays the centrality dependence of the V 2 -V 3 correlation (thick solid line) with those 
obtained for different q 2 event classes (markers). The correlation within a fixed centrality interval follows a path 
very different from the centrality dependence: the V 2 and V 3 are always anticorrelated with each other within a given 
centrality, whereas they are positively correlated as a function of centrality. Since the V 2 and V 3 are driven by the 
initial eccentricities, V2 oc 62 and V3 oc e 3 , one may expect similar anticorrelation between and £3- Indeed, a 
calculation based on a multiphase transport model [46] shows that such anticorrelations exist in the initial geometry 
and they are transferred into similar anticorrelations between V 2 and V 3 by the collective expansion [28]. 
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FIG. 6: (Color online) The correlation of v m between the 0.5 < pt < 2 GeV (x-axis) and 3 < pt < 4 GeV range (p-axis) for 
m = 2 (left panel) and m = 3 (right panel) in various centrality intervals. The data points are calculated in various q m intervals 
defined in Fig. 1 for each centrality, and they increase monotonically with increasing q m value. These data are overlaid with 
the centrality dependence without q m selection from Fig. 5. The thin solid straight lines represent a linear fit of the data in 
each centrality interval, and error bars represent the statistical uncertainties. 



FIG. 7: (Color online) The correlation of V 2 (x-axis) with V 3 (y- axis) both measured in 0.5 < pt < 2 GeV. The left panel 
shows the V2 and V3 values for fourteen 5% centrality intervals over the centrality range 0-70% without event-shape selection. 
The data points are connected to show the boomerang trend from central to peripheral collisions as indicated. The right panel 
shows the V2 and V3 values in the fifteen 52 intervals in seven centrality ranges (markers) with larger V2 value corresponding to 
larger 52 value; they are overlaid with the centrality dependence from the left panel. The error bars and shaded boxes represent 
the statistical and systematic uncertainties, respectively. 


In order to illustrate this anticorrelation more clearly, the V 2 -V 3 correlation data are replotted in Fig. 8 , separately 
for each centrality. The data are compared with the 62-63 correlations calculated via Eq. (2) from the MC Glauber 
model [34] and the MC-KLN model [47]. The MC-KLN model is based on the MC Glauber model, but takes into 
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FIG. 8: (Color online) The correlation of V 2 (*-axis) with V 3 (y- axis) in 0.5 < pr < 2 GeV for fifteen 52 selections in thirteen 
5% centrality intervals. The data are compared with the rescaled £ 2-63 correlation from MC Glauber and MC-KLN models in 
the same centrality interval. The data are also parameterized with a linear function [Eq. (15)], taking into account both the 
statistical and systematic uncertainties. The N par t dependence of the fit parameters is shown in the last two panels. The error 
bars and shaded bands represent the statistical and systematic uncertainties, respectively. 
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account gluon saturation effects in the initial geometry. One hundred million events were generated for each model 
and grouped into centrality intervals according to the impact parameter. The RMS e n value for each centrality interval 
is rescaled by a factor s n to match the inclusive v n value, which effectively is also the RMS value of v n [see Eq. (13)]: 


Sn — 



(14) 


The parameter s n changes with centrality but is assumed to be a constant within a given centrality interval. These 
constants are then used to rescale the 62-63 correlation to be compared with the V 2 -V 3 correlation in each centrality 
interval, as shown in Fig. 8 . In most centrality intervals the rescaled 63-63 correlation shows very good agreement with 
the V 2 -V 3 correlation seen in the data. However, significant deviations are observed in more central collisions (0-20% 
centrality range). Therefore, the V 2 -V 3 correlation data presented in this analysis can provide valuable constraints for 
further tuning of the initial geometry models. The V 2 -V 3 correlations in Fig. 8 are parameterized by a linear function: 

V 3 = kv 2 + v° , (15) 

where the intercept v° provides an estimate of the asymptotic V 3 value for events that have zero V 2 for each centrality. 
The fit parameters are summarized as a function of centrality (-/V part ) in the last two panels of Fig. 8 . 


D. 


V2—V4 and V3—V4 


correlations 



FIG. 9: (Color online) The correlation of V 2 (x-axis) with V 4 (y-axis) both measured in 0.5 < pt < 2 GeV. The left panel 
shows the V 2 and 14 values for thirteen 5% centrality intervals over the centrality range 0-65% without event-shape selection. 
The data points are connected to show the boomerang trend from central to peripheral collisions as indicated. The right panel 
shows the V 2 and V 4 values in different 52 intervals in seven centrality ranges (markers) with larger V 2 value corresponding to 
larger 52 value; they are overlaid with the centrality dependence from the left panel. The error bars and shaded boxes represent 
the statistical and systematic uncertainties, respectively. 


Figure 9(a) shows the correlation between V 2 and V 4 in 0.5 < pt < 2 GeV prior to the event-shape selection. 
The boomerang-like structure is less pronounced than that for the V 2 -V 3 correlation shown in Fig. 7(a). Figure 9(b) 
shows the V 2 -V 4 correlation for different (72 event classes (markers) overlaid with the centrality dependence taken from 
Fig. 9(a) (thick solid line). The correlation within a given centrality interval is broadly similar to the trend of the 
correlation without event-shape selection, but without any boomerang effect. Instead, the shape of the correlation 
exhibits a nonlinear rise for large V 2 values. 

To understand further the role of the linear and nonlinear contributions to V4, the V2-V4 correlation data in Fig. 9 
are shown again in Fig. 10, separately for each centrality. The data are compared with the 62-64 correlation rescaled 
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FIG. 10: (Color online) The correlation of Vi (a>axis) with va (y-axis) in 0.5 < pt < 2 GeV for fifteen q2 selections in thirteen 
5% centrality intervals. The data are compared with the rescaled £ 2-£4 correlation from MC Glauber and MC-KLN models in 
the same centrality interval. The data are also parameterized with Eq. (16), taking into account both statistical and systematic 
uncertainties. The N pal t dependence of the fit parameters is shown in the last two panels. The error bars and shaded bands 
represent the statistical and systematic uncertainties, respectively. 
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FIG. 11: (Color online) The centrality (N pa , T t) dependence of the V 4 in 0.5 < pr < 2 GeV and the associated linear and 
nonlinear components extracted from the fits in Fig. 10 and Eq. (17). They are compared with the linear and nonlinear 
component estimated from the previously published event-plane correlations [14] via Eq. (18). The error bars represent the 
statistical uncertainties, while the shaded bands or hashed bands represent the systematic uncertainties. 


according to Eq. (14). The rescaled e 2 -C 4 correlations fail to describe the data, suggesting that the linear component 
alone associated with £4 in Eq. (5) is not sufficient to explain the measured V 2 -V 4 correlation. 

To separate the linear and nonlinear components in the V 2 -V 4 correlation, the data are fitted to the following 
functional form: 


V 4 = \Jcl + (cic|) 2 . (16) 

This function is derived from Eq. (5), by ignoring the higher-order nonlinear terms (those in “...”) and a possible 
cross-term that is proportional to (cos 4 ( 4/2 — ^ 4 )). The fits, which are shown in Fig. 10, describe the data well for 
all centrality intervals. The excellent description of the data by the fits suggests that either the contributions from 
higher-order nonlinear terms and (cos4(’F 2 — ^ 4 )) are small, or the cross-term is, in effect, included in the nonlinear 
component of the fits. The centrality (iV part ) dependence of the fit parameters is shown in the last two panels of 
Fig. 10. 

The Co term from the fits can be used to decompose V4, without <72 selection, into linear and nonlinear terms for 
each centrality interval as: 


v\ = c 0 , = \Jv\- cl . (17) 

The results as a function of centrality are shown in Fig. 11 (open circles and squares). The linear term associated 
with £4 depends only weakly on centrality, and becomes the dominant part of V4 for lV par t > 150, or 0-30% centrality 
range. The nonlinear term increases as the collisions become more peripheral, and becomes the dominant part of V 4 
for lVp art <120. 

Since the contributions of higher-order nonlinear terms are small, as suggested by the fits discussed above, the 
linear and nonlinear contributions can also be estimated directly from the previously published event-plane correlator 











(cos 4($2 — $ 4 )) from ATLAS [14]: 
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NL.EP / a /\\ L,EP / 9 / NL,EP\ o /-i q\ 

v A = u 4 (cos4($2 - $4)) , V = yvi~(v, 4 ) 2 . ( 18 ) 

Results for this decomposition are shown in Fig. 11 (the hashed bands labeled EP), and they agree with the result 
obtained from the present analysis. 

Figure 12(a) shows the correlation between V 3 and u 4 in 0.5 < pt < 2 GeV prior to the event-shape selection. The 
data fall nearly on a single curve, reflecting the similar centrality dependence trends for V 3 and u 4 [11]. Figure 12(b) 
shows the V3—V4 correlation for different (73 event classes (colored symbols) overlaid with the centrality dependence 
taken from Fig. 12(a) (thick solid line). A slight anticorrelation between V 3 and v 4 is observed, which is consistent 
with the fact that V4 has a large nonlinear contribution from V 2 (Fig. 11), which in turn is anticorrelated with V 3 
(Fig. 7). 



FIG. 12: (Color online) The correlation of V 3 (i-axis) with V 4 (y- axis) both measured in 0.5 < pt < 2 GeV. The left panel 
shows the V3 and i > 4 values in thirteen 5% centrality intervals over the centrality range 0-65% without event-shape selection. 
The data points are connected to show the boomerang trend from central to peripheral collisions as indicated. The right panel 
shows the V3 and w 4 values for fourteen 53 selections (the two highest <73 intervals in Fig. 1 are combined) in several centrality 
ranges (markers) with larger 1)3 value corresponding to larger <73 value; they are overlaid with the centrality dependence from 
the left panel. The error bars and shaded boxes represent the statistical and systematic uncertainties, respectively. 


E. V2—V5 and V3—V5 correlations 

The analysis of V2-V5 and V3 — 1)5 correlations proceeds in the same manner as for the V2-V4 and V3-V4 correlations. 
A separation of the linear and nonlinear components of U 5 is made. 

Figure 13 shows the V 2 -V 5 correlation in 0.5 < pt < 2 GeV with <72 selection, separately for each centrality interval. 
The data are compared with the £ 2-65 correlations rescaled according to Eq. (14). The rescaled £ 2-65 correlations fail 
to describe the data in all centrality intervals, suggesting that the nonlinear contribution in Eq. ( 6 ) is important. To 
separate the linear and nonlinear component in the V2—V5 correlation, the data are fitted with the following function: 

= \Jcl + (civ 2 u 3 ) 2 , (19) 

where the higher-order nonlinear terms in Eq. ( 6 ) and a possible cross-term associated with (cos( 2 *F 2 + 34/3 — S'ks)) 
are dropped. For each centrality interval, Eq. (15) is used to fix the V 3 value for each value. The fits are shown in 
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FIG. 13: (Color online) The correlation of vi (tr-axis) with vs (y- axis) in 0.5 < pt < 2 GeV for fourteen <72 selections (the 
two highest 52 intervals in Fig. 1 are combined) in nine 5% centrality intervals. The data are compared with the rescaled 
£ 2-65 correlation from MC Glauber and MC-KLN models in the same centrality interval. The data are also parameterized with 
Eq. (19), taking into account both statistical and systematic uncertainties. The N part dependence of the fit parameters is shown 
in the last two panels. The error bars and shaded bands represent the statistical and systematic uncertainties, respectively. 


Fig. 13 and describe the data well for all centrality intervals. The centrality (lVp art ) dependence of the fit parameters 
is shown in the last two panels of Fig. 13. The Cq represents an estimate of the linear component of V 5 , and the 
nonlinear term is driven by ci, which has a value of ~ 1.5-2. 

Figure 14 shows the v 3 -v 3 correlations with q 3 selection in various centrality intervals. If Eq. (19) is a valid 
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FIG. 14: (Color online) The correlation of V 3 (x-axis) with V 5 (y- axis) in 0.5 < pt < 2 GeV for fourteen <73 selections (the 
two highest 53 intervals in Fig. 1 are combined) in nine 5% centrality intervals. The data are compared with the rescaled 
£ 3-65 correlation from MC Glauber and MC-KLN models in the same centrality interval. The data are also parameterized with 
Eq. (19), taking into account both statistical and systematic uncertainties. The N part dependence of the fit parameters is shown 
in the last two panels. The error bars and shaded bands represent the statistical and systematic uncertainties, respectively. 


decomposition of v$, then it should also describe these correlations. Figure 14 shows that this indeed is the case. 
The parameters extracted from a fit to Eq. (19) as shown in the last two panels of Fig. 14 are consistent with those 
obtained from V 2 -V 5 correlations. 

From the fit results in Figs. 13 and 14, the inclusive v$ values prior to event-shape selection are decomposed into 
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FIG. 15: (Color online) The centrality (7Vp ar t) dependence of the V 5 in 0.5 < pt < 2 GeV and the associated linear and 
nonlinear components extracted from the fits in Fig. 13, Fig. 14 and Eq. (20). They are compared with the linear and nonlinear 
component estimated from the previous published event-plane correlation [14] via Eq. (21). The error bars represent the 
statistical uncertainties, while the shaded bands or hashed bands represent the systematic uncertainties. 


linear and nonlinear terms for each centrality interval as: 


v 5 ~ c o> 



( 20 ) 


The results as a function of centrality are shown in the two panels of Fig. 15, corresponding to Fig. 13 and 14, 
respectively. Results for the two decompositions show consistent centrality dependence: the linear term associated 
with 65 dominates the v$ signal only in the most central collisions (iV part > 300 or 0-10% centrality). The nonlinear 
term increases as the collisions become more peripheral, and becomes the dominant part of v§ for !Vp art < 280. 

Similar to the case of V 2 -V 4 correlation, the linear and nonlinear contribution to v$ can also be estimated directly 
from the previously published event-plane correlator (cos( 2$2 + 3$3 — 54 > 5 )) from ATLAS [14]: 


NL,EP / /o/r 1 o/Fv r/R \\ L,EP / 2 / NL.EP\o 

v 5 = v 5 (cos( 2$2 + 3$3 - 5$ 5 )), v 5 ’ = \Jv £-{v 5 )- . 

Results for this decomposition are shown as solid curves in Fig. 15, and they agree well with the 
the present analysis. 


( 21 ) 

result obtained in 


F. Eccentricity-scaled v n 

One quantity often used to characterize the collective response of the medium to the initial geometry is the response 
coefficient k n defined in Eq. (4). Since the v n obtained from the two-particle correlation method effectively measure 
the RMS values of the event-by-event v n [43], a more appropriate quantity to characterize the collective response is 
the ratio of v n to the RMS eccentricity [22, 24]: v n j (e£). This quantity can be directly calculated for v 2 and V3 since 
they are mostly driven by £2 and £ 3 . But for higher-order flow harmonics, it is more appropriate to use the extracted 
linear component tA to make the ratios as it is more directly related to the e n . The v„ is taken as the eg term obtained 
from the two-component fits in Fig. 10 for n = 4 and Fig. 13 for n = 5. Figure 16 shows the centrality dependence of 
v n/ \/(£„) for n = 2 and 3 and u^/-\/(e^) for n = 4 and 5 (denoted by “linear” in figure legend), with e„ calculated 
in the MC Glauber model (left panel) and MC-KLN model (right panel). The higher-order flow harmonics show 
increasingly strong centrality dependence, which is consistent with the stronger viscous-damping effects as expected 
from hydrodynamic model calculations [16, 48, 49]. For comparison, the ratios are also shown for the full iq and V 5 
values without the linear and nonlinear decomposition, i.e., iq /\J (e|) (open diamonds) and I's/y^I) (open crosses); 
they show much weaker centrality dependence due to the dominance of nonlinear contributions to more peripheral 
collisions. 
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FIG. 16: (Color online) The eccentricity-scaled v n or the estimated linear component v„ obtained from two-component fits, 
V 21 yj (e|) (circles), V 31 yj (e§) (boxes), v\ / y/ (e|) (solid diamonds), v\ /yj (e§) (solid crosses), v±/ yj (t\) (open diamonds), and 
v$/y/ (eg) (open crosses). The eccentricities are calculated from the MC Glauber model (left panel) and the MC-KLN model 
(right panel). The error bars represent the statistical uncertainties, while the shaded bands or hashed bands represent the 
systematic uncertainties. 


VI. CONCLUSION 

Correlations between v m coefficients for m = 2 or 3 in different px ranges, and the correlation between v m and 
other flow harmonics v n for n = 2 to 5 in the same px range, are presented using 7 pb _1 of Pb+Pb collision data at 
y/s NN = 2.76 TeV collected in 2010 by the ATLAS experiment at the LHC. The v m -v n correlations are measured for 
events within a given narrow centrality interval using an event-shape selection method. Beside the centrality selection, 
this method makes a further classification of events according to their raw elliptic flow signal <72 or raw triangular 
flow signal q 3 in the forward rapidity range 3.3 < |?y| < 4.8. For each q m bin, the v m and v n coefficients are calculated 
at midrapidity |?/| < 2.5 using a two-particle correlation method, and hence the differential v m -v n correlation within 
each centrality interval can be obtained. 

The correlation of v m between two different px ranges shows a complex centrality dependence, but within a narrow 
centrality interval the correlation varies linearly with the event shape as determined by 52 or 53 . This linearity 
indicates that the viscous effects are controlled by the system size not by its overall shape. An anticorrelation is 
observed between V 2 and V 3 within a given centrality interval, and agrees qualitatively with similar anticorrelation 
between corresponding eccentricities £2 and £ 3 , indicating that these correlations are associated with initial geometry 
effects. 

The rq is found to increase strongly with V 2 , and U 5 is found to increase strongly with both v? and V 3 within a given 
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centrality interval. The trends and the strengths of V2-V4, V2-V5, and V3-V5 correlations disagree with corresponding 
e m -e n correlations predicted by MC Glauber and MC-KLN initial geometry models. Instead, these correlations are 
found to be consistent with a combination of a linear contribution to V 4 from e 4 and to V 5 from £ 5 , together with 
a nonlinear term that is a function of or of V 2 V 3 , as predicted by hydrodynamic models [23, 26]. The functional 
form of these nonlinear contributions is eclipsed in the overall centrality dependence, but has been directly exposed 
in the event-shape selected measurements reported here. A simple two-component fit is used to separate these two 
contributions in V 4 and V 5 . The extracted linear and nonlinear contributions are found to be consistent with those 
obtained from previously measured event-plane correlations. 

In order to quantify the response of the medium to the initial geometry, the extracted linear components of V4 
and V5, v\ and Vg, are scaled by the RMS eccentricity of corresponding order. The scaled quantities, v\/ (e|) and 
v\/\J (e§), show stronger centrality dependence than the similarly scaled quantities for elliptic flow and triangular 
flow, consistent with the stronger viscous-damping effects expected for higher-order harmonics. 
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